arXiv: 1503.03134vl [astro-ph.CO] 11 Mar 2015 


Galaxy Groups: A 2MASS Catalog 


R. Brent Tully, 

Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822, USA 


ABSTRACT 

A galaxy group catalog is built from the sample of the 2MASS Redshift Survey almost complete 
to K s = 11.75 over 91% of the sky. Constraints in the construction of the groups were provided 
by scaling relations determined by close examination of well defined groups with masses between 
10 11 and 10 15 M@. Group masses inferred from K s luminosities are statistically in agreement 
with masses calculated from application of the virial theorem. While groups have been identified 
over the full redshift range of the sample, the properties of the nearest and farthest groups are 
uncertain and subsequent analysis has only considered groups with velocities between 3,000 and 
10,000 km s -1 . The 24,044 galaxies in this range are identified with 13,607 entities, 3,461 of them 
with two or more members. A group mass function is constructed. The Sheth—Tormen formalism 
provides a good fit to the shape of the mass function for group masses above 6/i -1 x 10 12 Mq 
but the count normalization is poor. Summing all the mass associated with the galaxy groups 
between 3,000 and 10,000 km s -1 gives a density of collapsed matter as a fraction of the critical 
density of D coUapsed = 0.16. 
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mass function 


1. Introduction 

The definition of galaxy groups is fraught with 
ambiguity. Practitioners of simulated N-body uni¬ 
verses have a relatively easy time, given precise 
knowledge of the three-dimensional positions and 
velocities of large numbers of test particles. With 
observations of the real universe only the two pro¬ 
jected dimensions of galaxy positions are known 
with high accuracy. The third dimension, dis¬ 
tance, is usefully measured only for the nearest 
systems. Radial velocities are only crude discrim¬ 
inants, strongly modified by non-linear gravita¬ 
tional effects. Of proper motions we generally 
know nothing. Then there are the problems of 
small numbers and lost information. Most galax¬ 
ies live in groups that contain only a few members 
substantial enough to be recorded at a significant 
distance. As distance increases the recorded frac¬ 
tion of group members drops. 

The challenge to a cataloger is to capture the 
identities of groups consistently over wide ranges 


of both group mass and distance. Commonly, 
group constructio ns are built with a f riends-of- 
friends algorithm ( Huclira fe Gelleilll982lh The in¬ 
formation is pairwise separations in projection and 
radial velocity. The spatial separation scale chosen 
for the group construction is clearly dependent on 
distance since the density of recovered candidates 
falls with distance. Simultaneously, physically ap¬ 
propriate separation scales vary with group mass. 
Similarly, but even more sensitively, group veloc¬ 
ity differentials depend on mass. It is possible to 
incorporate the additional information of galaxy 
luminosities into the linkages. A recipe relating 
light to mass can be introduced. 

A catalog can be built under a set of assump¬ 
tions about the coupling of separations in space 
and velocity across large ranges in mass and dis¬ 
tance. Different assumptions may achieve inclu¬ 
sion of true group members and rejection of inter¬ 
lopers with varying degrees of success. Ultimately 
at issue is to what degree the properties of the cat¬ 
aloged groups are a correct representation of na- 
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ture as opposed to a biased representation imposed 
by the assumptions within the group-finding algo¬ 
rithm. There have been many attempts to build 
group catalogs si nce a n early qualitative effort by 
de Vaucouleursl ( 1975l f. Many i nvolve variations 


of _ fri ends-of-friends linkages dHuchra fe Gellerl 

1982: Geller fc Huchra 198. '4 Maia et al. 19891: 


rent group catalog are drawn from the 2MASS 
Redshift Survey (2MRS) sample that, aside from a 
narrow low Galactic latitude exclusion zone, is a l¬ 
most co mplete to K s = 11.75 ( Huchra et all 2 0l3) . 
Crook et al. ( 2007lf constructed a group catalog 


Ghnldn G^dj200^ lamella et al.|2002t Merchan fc 

20021 : lEke et al.l2no4ICrook et al.l2007t Lavau x fc Hudson skv. 


with a 2MR.S sam ple limited at K s = 11.25 and 
Lavaux fc Hudsonl ( 2011 1 built a 2MASS catalog 


ux limits that differ across regions of the 


2011: Temp el et all 2014lf. Others construct 


trees and/or use luminositi e s as a criterion for 


connectivity ( Mat,erne| 1978 1 _ C3miik|_ 1984; 


1987t iNolthenius fc White 1987 "[Nolthenius 


Tiilly 


1993; 


Yang et al 2*005^ 1 Weinmann et aH 20061: iMakarov Sz Karac^^ij§ey^ n 


20 llll . This new catalog is built with a methodol- 


The methodology that will be discussed in¬ 
volves a translation from luminosities to masses. 
Since the input catalog is flux limited, an adjust¬ 
ment is required to account for lost light with dis- 
galaxies in the 2MRS 11.75 sample are 


ogy th at most closely resembles that bvlYang et al 
(l2005lf . 

The present project is motivated by the propo¬ 
sition that we have a basic understanding of 
what constitutes a group so the search parame¬ 
ters should conform to this knowledge. Clearly 
there is a risk that the constraints will prejudice 
the outcome. However it can be argued that all ex¬ 
tant group catalogs have been prejudiced by their 
choices of constraints. Indeed in all extant group 
catalogs it is not hard to find proposed groups 
that are questionable. A test of the current cata¬ 
log will be whether it is a better survivor of this 
critique. 

A companion paper (Tully I2015TI . hereafter 
T15, presents the rationale for the group-finding 
algorithm that will be used. In that study, a mod¬ 
est number of very well studied groups are given 
attention. Those groups range from the Coma 
Cluster at the high mass end, at ~ 10 15 M 0 , to 
associations of dwarf galaxies with masses of a 
few times 10 11 A7q. Some of the groups are over¬ 
whelmingly dominated by early-type HI gas-poor 
systems and others are dominated by late-type 
gas-rich systems. In these well studied cases there 
are signatures of the extent of the (quasi) viri- 
alized halos and infall regions. It was possible 
to quantify scaling relations between dimensions, 
velocity dispersions, and masses. It is found that 
mass- to-light ratios for groups incre ase with group 
mass (iTullvl 20051 Yang et, akll2009l f. These group 
characteristics seem well enough established that 
they ought to be used as constraints in the con¬ 
struction of a group catalog. 

The galaxies that will be considered for the cur¬ 


given group assignments but the physical nature 
of the groups at expansion velocities above 10,000 
km s -1 are suspect. The lost light adjustment 
flairs to large values at large distances, begging 
the imposition of an upper limit to the practical 
catalog at 10,000 km s -1 in the CMB frame 
There is also a practical low velocity limit of 
3,000 km s -1 (CMB frame). Relative group pa¬ 
rameters are determined from redshifts (assuming 
H 0 = 100Y km s -1 Mpc -1 ) and peculiar veloc¬ 
ities add scatter to parameters at low observed 
velocities. Also, there are interesting options to 
a 2MRS sample locally. 2MASS samples misses 
low surface brightness galaxies. Such galaxies do 
not make an important contribution beyond 3,000 
km s^ 1 b ut are do minant in a nearby sample 
( Karachentsev et al.l 12013 ). At present, perhaps 
the best group catalog for the volume within 3,000 
km s -1 is the one by iMakarov fe Karachentsev) 


( 201 lh . The ambition here is to provide the best 
group catalog for the volume between 3,000 and 
10,000 km s -1 . All galaxies in 2MRS 11.75 are 
given an assignment reflecting their environment, 
be that to a group of two or more members or as 
a “group” of one. 

The outline of this paper is as follows. In the 
next section there will be a review of what can 
be claimed to be understood about the properties 
of galaxy groups, drawing heavily on the compan¬ 
ion paper T15. The third section will discuss the 
2MRS 11.75 base catalog. There will be an eval¬ 
uation of completeness issues, then a formulation 
of a group-finding algorithm consonant with the 
ideas presented in section 2. Section four hosts 
the catalog. In section five there is a discussion of 
the statistics of group properties. A group mass 
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function is constructed in section six. 


2. Properties of Known Galaxy Groups 

Dimension parameters that are accessible to N- 
body simulators, like the overdensity scale r 2 oo, a 
radius at 200 times the density required to close 
the universe with matter, or the dynamical gravi¬ 
tational radius r g 

r _ /-g 

E/>/ TfliTflj j T jj 

are not very useful to observers. How does one 
define a density surface about groups with only 
a few members? Likewise, the projected gravita¬ 
tional radius is unstable if statistics are poor and, 
in addition, requires a postulation of group mem¬ 
bership to be applied. 

It was argued in the companion paper T15 
that a group manifests a scale that is accessible 
to observations: the projected radius of second 
turnaround, i? 2 t (the nomenclature of T15 is fol¬ 
lowed, with projected dimensions in upper case 
R and those in 3D in lower case r). With ra¬ 
dial infall in a spherical system, galaxies initially 
plunge through the group and then continue out¬ 
ward until they stall at the second turnaround 
radius, rgt- Simplest models anticipate a cusp 
at ri t ( Bertschingeilll985tlShandarin fe Zeldovich 
1989lh With N-body collapse models the cusp 
is damp ed b ut a density discon tinuity feature re¬ 
mains ( Vogelsberger et al 112009 ). 

The discontinuity associated with Rat can be 
seen in judiciously selected groups. Here “judi¬ 
cious” means the chosen groups have been stud¬ 
ied individually in detail, both with imaging and 
spectroscopy, to depths sufficient to generate large 
samples. Efforts were made to minimize interloper 
problems. Details are given in T15. 

A second parameter readily measured for these 
well studied groups is the radial component of the 
velocity dispersion a p for galaxies within the ra¬ 
dius i? 2 t- In T15 it is demonstrated that there 
is a linear relationship between the two measured 
parameters 0 


Up/Rat = 491 ± 11 h km s 1 Mpc 


( 2 ) 


1 The scaling relations cited in T15 are derived with distance 

measures consistent with Ho = 75 km s _1 Mpc - 1 . In this 

paper, distances are converted from velocity units assuming 

Hq = 100 km s 1 Mpc -1 and the scaling relations are 


established from consideration of 13 groups rang¬ 
ing from the M31 halo to the Coma Cluster. 


Given a census of group members within Rat , 
the projected gravitational radius can be calcu¬ 
lated. Mass weighting, as explicit in Eq. [I] is not 
a good idea because a mass-weighted solution is 
dominated by a small number of objects and is un¬ 
stable. The information from dwarf test particles 
is ignored. Instead, the unweighted formulation is 
used 


R 


9 


N 2 


E,;> ? T /R 


ij 


(3) 


where there are N group candidates and pair pro¬ 
jections are counted only once. With this informa¬ 
tion, the group virial mass, M v , can be calculated 


M v = oj D r g /G = (onr/2G)a p R g (4) 


where statistically r g = (n/2)R g and ct^d = 
s/oLOp. The parameter a depends on the nature 
of the orbits and in T15 there is justification for 
the choice a = 2.5. 

With a dynamical measure of M v in hand (or 
M 12 = M v / 10 12 M q ), T15 demonstrated the cor¬ 
relation 


R 2t = 0.178 ± 0.005/i _ 2/3 M^2 3 Mpc (5) 


with 8 well studied groups ranging in mass from 
2 x 10 12 Mq to 2 x 10 15 Mq. An alternate, not 
independent, formulation is 

M v = 1.5 h- 1 x 10 6 cr 3 * . (6) 


How can this information be used to define 
groups from a redshift catalog? Given knowledge 
of group mass then the projected group scale Rat 
and velocity dispersion a p would be statistically 
known. Initially there is no knowledge of group 
mass but there is information from a proxy: K 
band luminosity. It is next necessary to evaluate 
the relation M v /Lk for groups over a range of 
masses. 


As a starting point, there is a known depen¬ 
dence of mass-to-light rati o on group mass at 
optical bands. iTullv ( 2005h found the depen¬ 
dence in blue light Lb oc M(?- 7 at masses above 
10 12 Mq. Group luminosities drop off rapidly be¬ 
low 10 12 Mq but this is a regime that will not 


modified accordingly. Following the standard convention, 
h = Hq/100. 
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be o f concern here. iM arinoni fe Hudsonl |2002) 
and xan den Bosch et, al.l (J2003) reached similar 
conclusions indirectly from the need to reconcile 
the galaxy luminosity function and the halo mass 
function. 

At K s band the Coma Cluster provides a robust 
mass-to-light calibration at high mass. T15 gave 
attention to a deeper survey than 2MRS 11.75 , the 


2MA SS Extended Source Catalog (jJarrett et al 


2000h with 367 group redshifts. At a distance 


given by the group mean CMB frame redshift of 
7331 km s -1 and Ho = lOO/i km s -1 Mpc -1 (as 
will be assumed consistently in this discussion) 
the Coma Cluster projected second turn-around 
radius is i? 2 t = 2.2h -1 Mpc, the projected ve¬ 
locity dispersion is 971 ±51 km s , the group 
virial mass is M v = 1.7h -1 x 10 15 M 0 , and the 
K s luminosity is Lk b = 1.4 h~ 2 x 1O 13 L 0 , whence 
M v /Lk s = 120 HMq/Lq. Here, Lk b has received 
adjustments for incompletion as will be discussed 
in the next section. 

At the other extreme of the mass range of inter¬ 
est, there is good information from the numerical 
action modelin g of galaxy orbits in the vicinity of 
nearby groups ( Shava fe Tullvll2013 1. A character¬ 
istic value for the dominant spirals in small groups 
is M/Lpc a = 4O/i.M 0 /L 0 . Elliptical galaxies have 
25% more mass per unit K s light. 

This difference in M/L over the interval 10 12 — 
1O 15 M 0 is consistent with the dependence found 
bv lMakarov fe Karachentsev! 1 2011 ) of M„ oc L 1 ^ 15 


wher e M p is the projected mass (IHeisler et al 
1985 ). These consistent constraints motivates 
the following formula for the translation from ob¬ 
served luminosity to expected group mass M^ xp 


M* xp = 43 x lO iU (Xi O ) 11& M 0 


(7) 


where L io is the K s luminosity in units of 1O 1O L 0 
and assuming Ho = 100 km s -1 Mpc -1 . This 
equation leads to M v / Lk b = 43 M 0 / Lq at Lk b = 
lO lo I/ 0 and 121 Mq/Lq at Lk b = 1O 13 L 0 . Below 
and above these luminosity limits M v / Lk b is held 
constant at 43 and 121 Mq/Lq respectively. The 
uncertainty in this transformation from luminos¬ 
ity to mass is estimated at 20% (68% probability) 
from comparisons between virial and luminosity 
masses discussed in Section [5j 


3. The 2MRS 11.75 Catalog 


The 2MASS Redshift Survey (2MRS) with red- 
shifts available for 98% of galaxies from the Two 
Micron All Sky Survey ( Jarrett et al.l 1200(1 12003 1 
brighter than K s = 11.75 outside a Milky Way 
exclusion zone of 9% of the sky is the last_ publi¬ 


cation and tribute to John Huchra (Huchra et al 


201 2b The redshift catalog used here contains 
43,065 entries. In the following, 2MRS total mag¬ 
nitudes corrected for extinction are used, resulting 
in a completion limit that is slightly boosted to 
K, = 11.5. 


3.1. The 2MRS Luminosity Function 

As with any flux limited sample, intrinsically 
faint galaxies are lost with increasing distance, so 
a first order of business is to determine a selection 
function that describes this loss of information. 
Others have considered this problem with related 
data sets. Here, a blind formulation is derived and 
then compared with the literature. 

It is to be appreciated that the 2MASS cata¬ 
log misses low surface brightness galaxies. The 
2MASS exposures were short, observing with 
modest sized telescopes, against high and variable 
auroral sky emission in the K s band. The photo¬ 
metric integrity of the survey is exceptional but 
Figure 1 provides a demonstration, in comparison 


with a deeper survey of a small region (jTull v et al 


19961 ). of the loss of low surface brightness flux. 
Very low surface brightness galaxies are lost en¬ 
tirely from the 2MASS catalog. 

Although Figure [T] looks scary, the loss of 
low surface brightness galaxies is not a sub¬ 
stantial problem for the present study because 
a lower redshift limit is set for the ‘practical’ 
group catalog that is being assembled (the regime 
3,000—10,000 km s -1 ). Low surface brightness 
galaxies are numerically dominant in a local 
volume limited sample to faint magnitudes but 
they make a small contribution to integrated 
flux. Here, a 2MASS K s band luminosity func¬ 
tion is constructed to a limiting magnitude of 
Mk b = —19+51og/i. The contribution from galax¬ 
ies fainter than this limit would augment the to¬ 
tal luminosity budget by 2-5% depending on the 
steepness of the faint-end luminosity function. Of 
slightly greater concern is the missing flux from 
low surface brightness galaxies above the rnagni- 
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Fig. 1.- Demonstration of lost flux from low sur¬ 
face brightness galaxies in the 2MASS catalog for 
galaxies in the nearby Ursa Major Cluster. On 
the ordinate, the magnitude difference is shown 
between the 2MASS K s magnitude and the mag- 
nit udemeasured in the more sensitive observations 
bv iTullv et al.l ( 1996lh The exponential disk cen¬ 
tral surface brightness of a galaxy is given on the 
abscissa. Symbols in red represent gala xies that 


appear in the Large Galaxy Atlas (j.Iarrett et al 


2003) which is more rigorous than the Extended 
Source Catalog ( Jarrett et al. 2000). The crosses 
identify lower limits to the magnitude differentials 
in cases of galaxies excluded from the 2MASS cat¬ 
alog. The 2MASS survey looses flux with a strong 
dependence on surface brightness. 


tude cutoff. As a rough estimate, as much as 20% 
of the light in stars in galaxies might be missing 
from our inventory. 


In defining the luminosity function it is nec¬ 
essary to contend simultaneously for the loss of 
faint galaxies with distance and the uneven con¬ 
tributions with distance due to clumping in the 
distribution of galaxies. In the regime where pe¬ 
culiar velocities are only a small perturbation on 
expansion velocities, it is possible to deal with the 
clumping problem by fitting the luminosity func¬ 
tion complete to a limit defined by the inner red- 
shift boundary of the shell. A parameterized func¬ 
tion can be fit to the observed luminosity distri¬ 
bution to the completeness limit for a sequence of 
shells over a domain relevant to the input catalog. 
If the parameters of the fits are stable over suc¬ 
cessive shells then it can be considered that the 
luminosity function is reasonably defined. 


A Schechter function ( Schechteilll976ll is fit fol¬ 
lowing standard conventions. The function for ab¬ 
solute magnitude Mx s is of the form N{Mx s ) = 


N K exp[—10~ OA ^ MKs ~ M ^]10~ OA< '°‘ K+1 ^ MKs ~ M ^' > 

( 8 ) 

and has three fitting parameters: a faint end power 
law slope parameter ax, a bright end exponential 
cutoff parameter and a sample normalization 
parameter N x. Only the first two parameters are 
important for the present discussion. 


It has to be entertained that galaxy luminosity 
functions may vary with environment. Detailed 
studies suggest that the variations are, in fact, 
quite modest (Tully [2011). Variations in a at r 
band range from — 1.35 for dense, gas-poor envi¬ 
ronments to —1.2 for less evolved gas-rich environ¬ 
ments. There is a coupling in the parameter fits: 
a brighter choice of M* leads to a more negative 
choice of a. There is tentative evidence at optical 
bands that there is a weak environmental depen¬ 
dence on luminosity functions but the more robust 
conclusion is that any dependence is at most small. 


Before giving attention to the full 2MRS catalog 
in redshift shells, three nearby clusters are consid¬ 
ered that completely sample to the Mx, = —19.0 
limit: the Virgo, Fornax, and Ursa Major clusters 
(here, luminosities have been transformed to be 
consistent with H 0 = 100). With the Virgo Clus¬ 
ter alone, a sample of 150 galaxies is fit with the 
Mx and ax parameters —23.7 and —0.97. Ex- 
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tending the sample to include 234 galaxies in the 
three clusters gives —23.4 and —0.83. In com¬ 
paring these numbers, recall the cautionary note 
above regarding parameter coupling. 

Turning to the full sample, consideration is 
given to fits in 15 separate shells at half magni¬ 
tude increments of M l ^ m , describing the faintest 
galaxies accessed in shells over the range —19 < 
M 1 ^ 71 < —26. The corresponding inner velocity 
limits range from 1259 km s -1 to 30,000 km s -1 . 
(With completion to K s = 11.5 then there is a 
complete representation of the luminosity function 
to Mk s = —19 + 51og/i for the volume within ex¬ 
pansion velocities of 1259 km s _1 .) The nearest 
shells contain ~ 500 galaxies. The most popu¬ 
lated shell at 10,000 km s -1 contains almost 8,000 
galaxies. 

The best constraints on the faint end slope ax 
come from fits to the nearest shells. However in 
these cases the bright end is poorly populated, cre¬ 
ating uncertainty in M^. The bright end expo¬ 
nential cutoff parameter is most confidently 
constrained in the mid range of shells which are 
well populated. In the five most distant shells be¬ 
yond 10,000 km s -1 , M l ^ m < so ax is effec¬ 
tively unconstrained and the fits are affected by 
rare high luminosity systems. 

Giving attention to the nearest shells there is 
the hint that the luminosity function is turning 
over (aji < — 1) but the evidence is not convinc¬ 
ing. This hint was already seen with the fits to 
the three nearby clusters. However the details of 
the faint end slope are not important here. At the 
near limit of 3,000 km s -1 of the ‘practical’ group 
catalog that will be constructed, = —20.9. 

For the ensuing fits ax — —1.0 is assumed and the 
only free parameter is M From the discussion 
earlier in this section it is to be appreciated that 
the faint end slope measurement for the 2MASS 
catalog is strongly biased shallow by the loss of 
low surface brightness galaxies. 

Figure [2] shows the run of fits to as a func¬ 
tion of the magnitude limit Mj£ m in shells, assum¬ 
ing ak = —1.0. In the four most distant shells, 
at velocities > 16,000 km s , there is a drift in 
M£ to higher luminosity. These shells only sam¬ 
ple the brightest galaxies in the exponential cutoff 
tail. There is a hint here that the Schechter de¬ 
scription is inadequate, that there is an excess of 
extremely luminous galaxies. Indeed, there is evi¬ 


dence from a more extensive sample derived from 
a larger volume that there are more bright galaxies 
than anticipated by the exponential cutoff of the 
Schechter function and that a better descr iption 


is provided by a double power law (Tempel et al 


2009il . That possibility is not of concern to the 


current discussion given an upper expansion veloc¬ 
ity cutoff of 10,000 km s -1 to the practical group 
catalog. 



Fig. 2.— Best fits for M* K in discrete redshift 
shells assuming ax = —1.0. Expansion velocities 
at the inner edge of shells are given in units of 100 
km s -1 along the bottom of the figure. Shells at 
velocities above 10,000 km s -1 only sample galax¬ 
ies brighter than M^. 


In summary of this section, the luminosity func¬ 
tion is adequately fit by the exponential cutoff pa¬ 
rameter M = —23.55 ± 0.05 + 51og/i as sumini 


ax = —1.0. By comparison, Cole et al. (2001 


cite —23.44 ± 0.03 with ax = 0.96 ± 0.05 from 
17,173 gala xies w ith redshifts from the 2dF survey, 
iBell et aD ( 2003 ) cite —23.33 with ax = —0.88 
brightward of —21 + 51og h and ax = —1.33 faint- 
ward, from 22,679 gal axies w ith redsh ifts from the 
Sloan survey, and Crook et al. ( 20071 ). referencing 
Huchra, cite —23.5 with ax = —1.02 from the 
preliminary 2MRS 11.25 sample. These earlier re¬ 
sults are consistent with the present results given 
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the issue of the lost low surface brightness systems 
and the coupling between the parameters and 
olk- 

3.2. The Selection Function 

Given a luminosity function, the number and 

flux from missing galaxies can be determined as 
M lrm 

increases in brightness with distance. Fig¬ 
ure [3] illustrates this lose with expansion velocity. 
Contributions faintward of Mk s = —19.0 + 51og h 
are ignored so there is complete coverage at 1259 
km s -1 . The correction factor is a multiplier to 
account for missing information. The number of 
missing galaxies explodes quickly but the lose of 
light increases much less rapidly. The luminosity 
correction factor CFi um is described by the poly¬ 
nomial expression 

CF, um = 1 + 1.68 x 10- 8 04 s - 1259) 2 

+2.2xl0” 17 (F L 5-1259) 4 +1.9xl0- 33 (V LS -1259) 8 

(9) 

where veloc i ties V t.s are in the Local Sheet frame 
f Tullv et al.lhoOSlb a variation on the more famil¬ 
iar but ambiguous Local Group frame. By 10,000 
km s , where M l £ m ~ M the correction factor 
for lost light is about a factor two. 

A procedure will be described in the next sec¬ 
tion that requires luminosities for the definition of 
galaxy groups. If CFi um > 2 the composition of 
groups become suspect with the proposed proce¬ 
dure. Consequently, although the full 2MRS 11.75 
catalog is given attention and groupings are pro¬ 
posed for the ensemble, group compositions are 
not considered reliable at velocities greater than 
10,000 km s _1 . At the other extreme, because 
of issues mentioned earlier associated with pecu¬ 
liar velocities and missing low surface brightness 
objects, the current catalog is not the best below 
3,000 km s _1 . The domain of the ‘practical’ cat¬ 
alog is 3,000—10,000 km s _1 . The properties of 
groups outside this regime are not considered re¬ 
liable. 


CF,„„= 1 + 1.68x 10»(V- 1259)' + a.2xlO-"(V-1259)'+ l.Bx 10-“(V- 1859)* 



Fig. 3.— Number and luminosity correction fac¬ 
tors as a function of velocity in the Local Sheet 
reference frame. Eq. [9] for CFi um is described by 
the solid line, overlying a dotted line that plots 
the numerically calculated run of this parameter. 


4. Construction of the Group Catalog 

Observable properties of established galaxy 
groups were discussed in section 2, including sec¬ 
ond turnaround dimensions i? 2 t distinguishable 
in well studied cases, the line-of-sight velocity 
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dispersion a v in clean cases, and the associated 
luminosities. Scaling laws were established that 
permit the inference of halo properties — mass, 
velocity dispersion, and radius - from observed 
luminosities. Section 3 provided a description of 
adjustments that have to be made to luminosities 
to account for lost contributions as a function of 
distance. 

The conceptual outline of the group-finding al¬ 
gorithm is as follows. (1) Start with the intrinsi¬ 
cally most luminous galaxy in the sample after ad¬ 
justment with the correction factor. (2) Assume 
a group mass-to-light ratio appropriate for that 
intrinsic luminosity using Eq. 0 and calculate the 
halo expectation parameters i? 2 t (Eq. 0 and a v 
(Eq. 0). (3) Cycle through the sample to search 
for galaxies that he within the R- 2 t radius of the 
primary system and within 2a p of its velocity. (4) 
After this first cycle, sum the luminosities of as¬ 
sociated galaxies and determine their luminosity 
weighted projected centroid and unweighted veloc¬ 
ity mean. (5) Calculate the halo expectation pa¬ 
rameters i? 2 t and &p for this enlarged entity. Re¬ 
peat cycles until there are no new links. (6) Go to 
the next intrinsically most luminous galaxy among 
the unlinked cases and repeat procedures (2) — (5). 
(7) Continue to successively fainter galaxies until 
there are no more galaxies to consider. 

After the initial construction of tentative 
groups, it is found that there can be occasional 
overlaps between close neighbors. Hence another 
cycle is initiated. (8) Beginning with the most 
populated candidate group, cycle through the 
other groups looking for overlaps in both i? 2 t pro¬ 
jected dimensions and velocity dispersions (the 
larger of the quadrature addition of the 2 a p val¬ 
ues and the 3cr p value of the larger entity). (9) 
Recalculate halo properties for any enlarged can¬ 
didate group and recycle. Repeat until no new 
additions. (10) Consider next most populous can¬ 
didate group and successively smaller candidates 
until the entire catalog has been explored. The 
final affiliations of galaxies will be called ” nests”. 

An illustration of groups found in a particular 
region is given in Figure 0 The Perseus-Pisces 
filament that is shown is the densest region of ma¬ 
jor filaments in the volume that has been explored. 
Histograms of the velocities of candidate members 
in the four largest nests are seen in Figure 0 

Among details, one is a cutoff in the galaxies 
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Fig. 4.— The Perseus-Pisces filament. Major 
components are given distinct colors. Dotted cir¬ 
cles indicate R 2 1 , second turnaround radii for the 
major halos. Points in black identify all other 
2MRS K < 11.75 galaxies with 3500 <V< 6500 
km s -1 . 
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Fig. 5.— Velocity histograms for the four most 
important components of the Perseus-Pisces fila¬ 
ment. 


that are considered at an expansion velocity of 
24,000 km s -1 . The idea is to identify groups 
in some manner to 20,000 km s -1 with mini¬ 
mal high velocity edge effect. As mentioned al¬ 
ready, though, the luminosity correction factor be¬ 
comes large at high redshifts, reaching 40 at 20,000 
km s -1 . The group catalog extends to 20,000 
km s _1 but it would be unwise to make much 
of group characteristics at velocities greater than 
10,000 km s” 1 . 

Another detail is that the north and south 
Galactic hemispheres were evaluated separately 
since no real group except the most local crosses 
the zone of obscuration. In the final group cat¬ 
alog, nests are given 6-figure identification num¬ 
bers with the initial digit for those in the north 
the number 1 and the initial digit for those in the 
south the number 2. The 2MRS limit of |6| = 5° 
fortuitously allows inclusion of two of the largest 
nearby clusters, Norma and Ophiuchus, as shown 
in Figure [6] The excluded low Galactic latitude 
zone occupies 9% of the sky. As an aside, a system¬ 
atic X-ray search for clusters has revealed 3 clus¬ 
ters within 10,000 km s _1 in the 2MRS exclusion 
zone: CIZAJ0450.0+4501 (14 = 6655 km s” 1 ) and 


CIZAJ0603.8+2939 (14 = 8994 km s" 1 ) in prox¬ 
imity to the Perseus Cluster and CIZAJ1324.7- 
5736 (14 = 5696 km s -1 ) near the Norma Cluster 
( Ebeling et al. 2002lh 

A third, more significant detail follows from 
close inspection of the nests in plots of projected 
positions and histograms of velocities. There 
sometimes were apparent gaps that suggest the 
linkages were too severe. These gaps were usu¬ 
ally more convincing in the spatial information 
rather than in velocities. In only six instances, 
nests were split. In these cases the smaller part is 
given an identification with the number 2 as the 
second digit. The digits 3 to 6 are the same as 
the identification number of the larger part. The 
Abell 1367 Cluster seen in Figure [7] is an exam¬ 
ple. The main cluster in red is nest 100005 while 
the secondary cluster in cyan is nest 120005 in the 
group catalog. It can be expected in such a case 
that the two entities are on paths to soon merge. 
Nests are split where the evidence is compelling. 
With several other nests there are hints of segre¬ 
gation but there is sufficient ambiguity that they 
were left alone. For example, in the case of the 
Centaurus Cluster shown in Figure [5] there ap¬ 
pear to be two kinematic components but they 
are not distinctly separated and the spatial over¬ 
lap is severe . The two components in this wel l 
known case ( Lucev et al.1 19861: Stein et al. 1997 1 
appear to have similar distances. Here and com¬ 
monly with the halos under consideration it can 
be expected that merging and cluster growth is 
an ongoing process. 


A fourth detail results from the observation 
that the algorithm that is employed fails to pick up 
extreme velocity outliers. With the large nests it 
was not uncommon to find velocity outliers reach¬ 
ing ~ 3.5cr projected very near to the center of 
the nest. Their central locations and the paucity 
of galaxies in the vicinity with similar velocities 
makes the likelihood of their membership very 
high. These galaxies are added to the associated 
nests and the nest parameters are recalculated. 

Table Q] gives the gross statistics of the group 
catalog, both overall and in the restricted inter¬ 
val 3,000 to 10,000 km s _1 . In the velocity in¬ 
terval 3,000 km s -1 to 10,000 km s -1 there are 
24,044 galaxies to consider in the 2MRS 11.75 
sample. These are assembled into 13,605 groups, 
with 3,461 groups of at least 2 members and 10,144 
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Nest S200002 PGC 57612 Norma 



Nest N100009 PGC 59827 Ophiuehus 
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Fig. 6.— Norma Cluster shown in the top panel 
and Ophiuehus Cluster shown in the bottom panel 
are two major clusters at low Galactic latitude. 
The plane of the Milky Way lies at SGL=180 and 
the grey lines conform to the boundary of the 
2MRS redshift survey. The PGC number iden¬ 
tifies the dominant galaxy in the cluster. 
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Fig. 7. The region around the Abell 1367 Clus¬ 
ter. The main cluster in red and an apparently 
distinct entity in cyan have strongly overlapping 
velocities. 


singles. The Perseus Cluster contains the most 
galaxies from the 2MRS 11.75 sample in this re¬ 
stricted velocity range, with 180. Figure [3] teaches 
us that number counts are strongly dependent 
on distance. Masses from adjusted luminosities 
and velocity dispersions are more stable with dis¬ 
tance. Table [2] lists the top 8 clusters in the 
3—10 thousand velocity interval: on the left, 
by masses derived from luminosities, and on the 
right, by velocity dispersions. For comparison, the 
Virgo Cluster, the most important cluster within 
3,000 km s -1 is included. The Virgo Cluster is a 
comfortable peer of these prominent groups. 

4.1. Local Densities 

Knowledge of the selection function and of the 
group characteristics allows for the determination 
of the smoothed luminous density at the location 
of every galaxy in the sample. Values for this 
parameter have been calculated on a l/i -1 Mpc 
grid and each galaxy is given the value at the 
nearest grid location. The luminosity densities at 
each grid point are given by the sum of contri- 
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Table 1 
Summary 


Sample 

ff All Galaxies 

# Groups to 20k 

> 2 

# 3k-10k 

# Groups 

> 2 

North 

21,995 

13,090 

3,178 

12,153 

6,982 

1,781 

South 

21,043 

12,376 

3,010 

11,891 

6,624 

1,680 

Total 

43,038 

25,475 

6,188 

24,044 

13,606 

3,461 


Table 2 

Most Prominent Clusters 

Group 

Name 

Mass h- 1 ! 0 15 

M q CF 

Group 

Name 

(7p km s 1 

100004 

A2199 

1.8 

2.3 

100009 

Ophiuchus 

976 

100001 

Coma 

1.7 

1.6 

200001 

Perseus 

962 

100009 

Ophiuchus 

1.6 

2.3 

200002 

Norma 

957 

200002 

Norma 

1.3 

1.3 

100001 

Coma 

886 

200001 

Perseus 

1.2 

1.2 

100003 

Centaurus 

822 

100002 

Virgo 

1.1 

1.0 

200017 

A539 

754 

200016 

A2634 

0.8 

2.3 

100004 

A2199 

740 

200017 

A539 

0.7 

1.9 

100002 

Virgo 

717 

100005 

A1367 

0.7 

1.5 

100005 

A1367 

707 


tuitions from members of the sample after Gaus¬ 
sian smoothing. Luminosities are adjusted follow¬ 
ing the recipe illustrated in Fig. [3] The Gaussian 
smoothing scale is 1 h~ l x CF ^ Mpc so the com¬ 
pensation for missing light is spread across more 
grid points at larger distances while the contribu¬ 
tion from a single galaxy remains constant at the 
central position of the galaxy. In order to elimi¬ 
nate the finger-of-god elongation of structures in 
redshift space, yet realistically replicate the depth 
of clusters, galaxies in the identified groups with 
velocity V ga i are given redshift space positions in 
the line-of-sight based on the mean group velocity 
V gp plus 0.1 (Vgai — V gp ). Finally to note: since 
the luminosity adjustment becomes very large at 
large redshifts, the luminosity density computa¬ 
tion is restricted to a box ±200ft -1 Mpc on a side. 
Even within this box, density values toward the 
edges are very uncertain. 

Maps of the density distribution are shown in 
three orthogonal slices in Figures [HI OH and HU 
The presentations are in supergalactic coordinates 
and each slice has a thickness of 4000 km s -1 . 


The slices are chosen to include most of the main 
structural features in the volume being studied. 
See the captions for details. 

4.2. Tabulated Information 

Three on-line tables are provided that present 
the group catalog. Table [3] gives a summary of 
each group, one line per group. Table [4] identifies 
the individual members of each group. Table [5] 
combines these pieces of information: there is one 
line per galaxy containing information about the 
galaxy and about the ensemble properties of its 
group. In each case, the first 10 rows of large ta¬ 
bles are shown, with the full tables made avail¬ 
able on-line. Table [5] is also available (separated 
into Galactic north and south halves) at the Ex- 
tragalactic Distance Database!! More detailed in¬ 
formation about the tables follows. 

Table [3' Summary of Group Properties. Column 
(1) Group (nest) identification number; lxxxxxx 


2 http://edd.ifa.hawaii.edu; see catalogs ‘2MRS North & 
South Groups’ 
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SGL 


Fig. 8.- The region around the Centaurus Clus¬ 
ter. Two relatively distinct kinematic components 
strongly overlap in projection. The blue and red 
components separated at 4300 km s^ 1 have been 
called Cen30 and Cen45 in the literature. 



Fig. 9.— Slice through smoothed 2MRS lu¬ 
minosity density cube: SGZ interval -2000 to 
+2000 km s -1 . Contours are at logarithmic inter¬ 
vals with a low density cutoff at 10 9 h 2 M 0 Mpc~ 3 . 
Galactic obscuration band at SGY ~ 0. The most 
prominent feature is the Shapley Concentration at 
SGX = -13,000 km s" 1 , SGY = 8,000 km s" 1 . 
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Fig. 10.— Slice through smoothed 2MRS lu¬ 
minosity density cube: SGY interval +5000 to 
+9000 km s _1 . Contours as in Fig. [9] This inter¬ 
val contains the Great Wall with the Coma Cluster 
near the center, the Hercules complex above center 
and slightly to the left, and the Shapley Concen¬ 
tration at SGX = —13,000 km s -1 , SGZ = 0. 
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Fig. 11. Slice through smoothed 2MRS lu¬ 
minosity density cube: SGX interval -4000 to 
0 km s -1 . Contours as in Fig. [9] Galactic ob¬ 
scuration is tilted by 6° from the SGY = 0 
axis. The Great Wall rises vertically at SGY = 
6, 000 km s^ 1 reaching the Hercules complex at 
SGZ + 8,000 km s -1 . Directly above our position 
at the origin is the Local Void, blending toward 
the upper right into the Hercules Void. 
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if in north galactic hemisphere and 2xxxxxx if in 
south. If the identification is 12xxxxx or 22xxxxx 
the group is a split from an adjacent larger group. 
(2) Number of members from the 2MASS 11.75 
catalog. (3) Principal Galaxies Catalog (PGC) 
name for brightest member in group. (4—5) Lumi¬ 
nosity weighted supergalactic longitude and lati¬ 
tude of group, in degrees. (6) Logarithm of K s 
band luminosity summed over all 2MASS 11.75 
members and adjusted by the correction factor for 
missing light, units of solar luminosity. (7) Un¬ 
weighted average group velocity, in the cosmic mi¬ 
crowave background frame adj usted by a cos mo- 
logical model as described in lTullv et, al. ( 2013 ). in 


km s -1 . (8) Distance modulus from group veloc¬ 
ity assuming Ho = 100 km s -1 Mpc -1 , in magni¬ 
tudes. (9) Line of sight bi-weight velocity disper¬ 
sion of group members; requires at least 5 group 
members, in km s -1 . (10) Projected radius of sec¬ 
ond turnaround calculated from group luminosity 
and assumed conversion to mass, at group distance 
given by velocity and Ho = 100 km s -1 Mpc -1 , in 
Mpc. (11) Line of sight velocity dispersion antici¬ 
pated from group luminosity and assumed conver¬ 
sion to mass, in km s — . (12) Group mass calcu¬ 
lated from luminosity, adjusted by correction fac¬ 
tor, and assumed light to mass conversion factor. 
(13) Conversion factor for missing light as a func¬ 
tion of systemic velocity. 

Table Group Members. Column (1) Group 
(nest) identification number; lxxxxxx if in north 
galactic hemisphere and 2xxxxxx if in south. If the 
identification is 12xxxxx or 22xxxxx the group is 
a split from an adjacent larger group. (2) Princi¬ 
pal Galaxies Catalog (PGC) name of group mem¬ 
ber. (3—4) Supergalactic longitude and latitude 
of galaxy, in degrees. ( 5) Ty pe i n the numeric 
code of Ide Vaucouleurs et al.l ( 199ll l. (6) 2MASS 
K s total magnitude, extinction corrected, in mag¬ 
nitudes. (7) 2MASS J — H total extinction cor¬ 
rected color, in magnitudes, (8) 2MASS J — K to¬ 
tal extinction corrected color, in magnitudes. (9) 
Logarithm of intrinsic K s luminosity, extinction 
corrected, at group distance given by velocity and 
Ho = 100 km s -1 Mpc -1 , in units of solar lumi¬ 
nosity. (10) Velocity of galaxy in the cosmic mi¬ 
crowave background frame a djusted by a cosmo ¬ 


logical model as described in iTullv et al 
in km s -1 . 


(2013|), 


Table El- All Galaxies with Group TagsE Col¬ 
umn (1) Principal Galaxies Catalog (PGC) name 
of galaxy. (2—3) Galactic longitude and latitude, 
in degrees. (4—5) Supergalactic longitude and lat¬ 
itude in degrees. (6) Type i n the numeric code of 
Ide Vaucouleurs et al.l ( 199ll ). (7) Heliocentric ve¬ 
locity, in km s -1 . (8) Velocity in_the L ocal S heet 


reference frame defined bv ITullv et al. ( 2008 ) 


km ° 1 


s \ (9) Velocity of galaxy in the cosmic mi¬ 
crowave background frame a djusted by a cosmo¬ 


logical model as described in Tully et al 


(1201311 . 

in km s -1 . (10) 2MASS J — H total extinc¬ 

tion corrected color, in magnitudes, (11) 2MASS 
J — K s total extinction corrected color, in mag¬ 
nitudes. (12) 2MASS K s total magnitude, ex¬ 
tinction corrected, in magnitudes. (13) Loga¬ 
rithm of intrinsic K s luminosity, extinction cor¬ 
rected, at group distance given by velocity and 
Ho = 100 km s -1 Mpc -1 , in units of solar lumi¬ 
nosity. (14) Logarithm of luminosity density at 
nearest location on a l/i -1 Mpc grid, where con¬ 
tributions come from all adjacent 2MASS 11.75 
galaxies after smoothing with a Gaussian of l/i -1 
Mpc times an adjustment dependent on the lumi¬ 
nosity correction factor as described in the text, 
in units of solar luminosity per Mpc 3 . (15) Group 
(nest) identification number; lxxxxxx if in north 
galactic hemisphere and 2xxxxxx if in south. If the 
identification is 12xxxxx or 22xxxxx the group is 
a split from an adjacent larger group. (16) Num¬ 
ber of members from the 2MASS 11.75 catalog. 
(17) Principal Galaxies Catalog (PGC) name for 
brightest member in group. (18—19) Distance and 
distance modulus from group velocity assuming 
Ho = 100 km s -1 Mpc -1 , in Mpc and magni¬ 
tudes respectively. (20—21) Luminosity weighted 
supergalactic longitude and latitude of group, in 
degrees. (22) Logarithm of K s band luminosity 
summed over all 2MASS 11.75 members and ad¬ 
justed by the correction factor for missing light, 
units of solar luminosity. (23) Conversion factor 
for missing light as a function of systemic velocity. 
(24) Line of sight velocity dispersion anticipated 
from group luminosity and assumed conversion to 
mass, in km s -1 . (25) Projected radius of sec¬ 
ond turnaround calculated from group luminosity 
and assumed conversion to mass, at group distance 
given by velocity and Hq = 100 km s -1 Mpc -1 , 


3 Available with updates at http://edd.ifa.hawaii.edu 
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in Mpc. (26) Unweighted average group velocity, 
in the cosmic microwave background frame ad¬ 
j usted by a cosm ological model as described in 


i u11v et aiJ ihnah in km s 1 . (27) Bi-weight 


group ([Beers et al. lH) velocity in the cosmic 
microwave background frame adjusted by the cos¬ 
mological model, in km s -1 . (28) Uncertainty in 
bi-weight group velocity, in km s -1 . (29) Line of 
sight bi-weight velocity dispersion of group mem¬ 
bers; requires at least 5 group members, in km s . 
(30) Line of sight velocity dispersion of group 
members, in km s — . Null if group of one. (31) Bi¬ 
weight projected gravitational radius Rij , in Mpc. 
(32) Uncertainty in bi-weight projected gravita¬ 
tional radius, in Mpc. (33) Group mass from 
virial theorem with bi-weight dispersion and ra¬ 
dius parameters, in units of 10 12 M 0 . (34) Group 
mass based on adjusted intrinsic luminosity and 
mass to light prescription, i n units of 10 12 M 0 . 
(35—36) ICrook et all ( 2007h high and l ow den ¬ 
sity gr oup identifications. (37) lLavaux fe Hudson! 
( 20111 1 group identification. (38—40) Supergalac- 
tic X,Y,Z coordinates from group velocities assum¬ 
ing Ho = 100 km s -1 Mpc -1 , with small adjust¬ 
ments described in section 14.11 to create roughly 
spherical clusters, in Mpc. 


5. Group Properties 


It can be asked if observed velocity dispersion, 
spatial scale, and inferred virial masses track the 
assumptions based on luminosity. Correlations are 
expected but there could be offsets. There will 
be evident uncertainties with small groups. In 
best cases with order 100 candidate members, pro¬ 
jected velocities and positions reasonably repre¬ 
sent the three-dimensional distribution but with 
only a few members, in the extreme only two, 
knowledge of only one of three velocity compo¬ 
nents and two of three spatial components results 
in large uncertainties in group properties. Un¬ 
certainties are compounded by interlopers. Then, 
more fundamental than the observational consid¬ 
erations, it is appreciated that the structures may 
often, even usually, stray from dynamical equilib¬ 
rium. 


While ultimately uncertainty will rule, there 
are better and worse ways to evaluate group pa- 
rame ters. Approaches have been discussed by 
Beers et al. ( 1990 1. They evaluate methods for 


what they call resistance, robustness, and effi¬ 
ciency. A resistant method is minimally affected 
by outliers. A robust method works with diverse 
population characteristics. An efficient method 
does the best that can be done with poor statis¬ 
tics. Following tests with samples from 5 to 200 in 
size, Beers et al. identify a preference for bi-weight 
location and scale estimators. Their recipes are 
followed here to determine group central veloci¬ 
ties and line-of-sight dispersions and measures of 
the projected separations Rij of group members. 

Figure [12] provides a comparison between group 
mass estimates derived via the virial theorem from 
positions and velocities versus group mass es¬ 
timates following from the integrated K s band 
light. The dotted line gives the equality relation¬ 
ship. It is seen that there is excellent agreement 
with groups containing at least 30 2MRS galaxies, 
identified by the brown and red points, and still 
good agreement with memberships as small as 10, 
points in orange and green. The correlation falls 
apart with fewer group members, already with 5 
to 9 members represented in the plot in cyan, but 
especially with groups of 2 to 4 represented in blue. 



Group Luminosity Mass [h 1 M 0 ] 

Fig. 12.— Comparison of group mass estimates 
based on luminosity vs. estimates based on veloc¬ 
ity dispersions and spatial separations. The color 
code identifies the number of members. 
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Table 3 

Nest Properties (1st 10 of 25474) 


Nest 

Mem 

PGC1 

SGL 

SGB 

log L b k 

VLd 

DM 

(T V 

R-2t 

<7p 

Mass 

CF 

100001 

136 

44715 

89.6226 

8.1461 

13.15 

7331 

34.33 

886 

2.393 

881 

0.171E+16 

1.65 

100002 

197 

41220 

103.0008 

-2.3248 

12.69 

1491 

30.87 

670 

1.617 

596 

0.529E+15 

1.00 

100003 

113 

43296 

156.2336 

-11.5868 

12.75 

3873 

32.94 

822 

1.708 

629 

0.623E+15 

1.12 

100004 

81 

58265 

71.5103 

49.7851 

13.16 

9424 

34.87 

740 

2.418 

891 

0.177E+16 

2.26 

100005 

61 

36487 

92.0255 

-10.4950 

12.78 

6987 

34.22 

707 

1.753 

646 

0.673E+15 

1.58 

100006 

85 

31478 

139.4478 

-37.6063 

12.48 

4099 

33.06 

648 

1.347 

496 

0.305E+15 

1.14 

100007 

86 

56962 

108.5182 

49.0878 

13.52 

11603 

35.32 

1261 

3.168 

1167 

0.398E+16 

3.30 

100008 

65 

39600 

67.2953 

3.2390 

11.94 

1054 

30.11 

209 

0.837 

308 

0.733E+14 

1.00 

100009 

66 

59827 

173.2412 

43.4150 

13.13 

9112 

34.80 

976 

2.352 

866 

0.163E+16 

2.15 

100010 

55 

47202 

149.1678 

-1.3439 

13.93 

15265 

35.92 

1002 

4.363 

1607 

0.104E+17 

7.95 


Table 4 

Nest Members (1st 10 of 43038) 


Nest 

PGC 

SGL 

SGB 

Ty 

Ktot 

J-H 

J-K 

log L k 

*mod 

100001 

44715 

89.6369 

8.3951 

-4.3 

8.40 

0.69 

0.99 

11.69 

6834 

100001 

44628 

89.6262 

8.2751 

-3.6 

8.86 

0.82 

0.98 

11.50 

7582 

100001 

43895 

90.2180 

6.9260 

-2.3 

9.24 

0.47 

0.75 

11.35 

8837 

100001 

44298 

89.9643 

7.6957 

-3.9 

9.19 

0.75 

0.88 

11.37 

7763 

100001 

44323 

89.0106 

7.9544 

-4.1 

9.22 

0.81 

1.03 

11.36 

7164 

100001 

44840 

89.8613 

8.5241 

4.0 

9.84 

0.67 

0.98 

11.11 

8420 

100001 

44938 

90.0745 

8.6928 

-3.1 

9.89 

0.67 

0.94 

11.09 

8300 

100001 

43981 

89.9140 

7.1398 

-2.6 

9.89 

0.70 

1.02 

11.09 

8288 

100001 

44737 

89.4237 

8.4824 

-2.1 

10.14 

0.68 

0.94 

10.99 

8998 

100001 

44114 

89.6574 

7.4958 

-3.0 

9.71 

0.68 

0.98 

11.16 

7323 


Table 5 

Combined Catalog (1st 10 of 43038) 


PGC 

Glon 

Glat 

SGL 

SGB 

Ty 

Vfiel 

V LS 

vmod 

J-H 

J-K 

Ktot 

log L k 

log PK 

2 

113.9554 

-14.6992 

341.6442 

20.7389 

3.1 

5017 

5309 

4799 

0.72 

1.04 

9.50 

10.91 

10.41 

5 

110.6206 

-28.9043 

326.1772 

19.7805 

-0.6 

10445 

10717 

10310 

0.70 

1.01 

10.49 

11.12 

10.34 

12 

90.1920 

-65.9300 

286.4249 

11.3510 

1.3 

6546 

6683 

6279 

0.72 

0.99 

11.11 

10.46 

9.36 

14 

101.7853 

-52.4728 

300.9143 

15.3673 

-4.9 

11602 

11801 

11571 

0.76 

1.03 

10.60 

11.22 

10.57 

16 

91.6006 

-64.8655 

287.6120 

11.7030 

-0.1 

5664 

5806 

5434 

0.63 

0.82 

11.52 

10.18 

9.78 

18 

113.9223 

-15.0086 

341.3140 

20.7012 

2.0 
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Of course, the virial mass estimates are suspect 
with small numbers. The systematic trend toward 
small masses can be understood. The group search 
algorithm that sets spatial and velocity limits fa¬ 
vors inclusion of galaxies with close projected sep¬ 
arations and radial velocities and disfavors objects 
with large projected separations and radial veloci¬ 
ties that would make the cut with full 3D informa¬ 
tion. The scaling laws assume that the accessible 
phase space information reasonably represents the 
group which is statistically valid if the member¬ 
ship is large. 

The separate velocity dispersion and spatial 
scale components are given attention in Figure ITdl 
Observed dispersions are compared with disper¬ 
sions anticipated from luminosities in the top 
panel and observed gravitational radii are com¬ 
pared with second turnaround radii in the bottom 
panel. Color codings tied to group membership 
are the same as in Fig. [12] except groups with 
less than 5 members are omitted. The absolute 
luminosity cutoff Mk s = —19 + 51og/i translates 
to lower cutoffs in the dispersions and dimensions 
derived from luminosities. The dashed lines in the 
two panels represent best fits to the comparisons 
involving groups with at least 30 members. In the 
case of velocity dispersions the difference between 
measures is marginal. In the case of dimensions, 
there is no expectation that the two measures 
would be the same. Statistically, i? 2t is 27% ± 3% 
larger than the gravitational radius. 

There are two conclusions from the discussion 
of this section. The first is that the absolute scale 
of the masses derived from luminosities is consis¬ 
tent with a scale from kinematic input. The sec¬ 
ond is that K s luminosities provide masses with 
smaller uncertainties than kinematic masses. Ad¬ 
mittedly there could be hidden systematics in the 
luminosity masses that remain to be revealed. 
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6. The Group Mass Function 

The availability of the group catalog presents 
the opportunity to construct a well formulated 
halo mass function. The analysis is restricted 
to groups with mean velocities in the range 
3,000—10,000 km s -1 in the CMB frame. The 
definition of groups includes groups of one, so the 
concepts of groups and halos are synonymous to 
the degree that halos are filled with enough stars 


Fig. 13.— Top: bi-weight radial velocity disper¬ 
sion vs. velocity dispersion inferred from K s lumi¬ 
nosity. Bottom: bi-weight projected gravitational 
radius vs. i? 2 t- The color code identifies the num¬ 
ber of members. Dotted curves correspond to 1:1 
dependencies. Dashed curves are the best fits to 
groups with at least 30 members. 
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to be visible. Within the completion limits of the 
2MRS 11.75 survey, the group catalog contains all 
halos in the volume. 

The mass function derived with this sample is 
presented in Figure [l4j Masses are determined 
from K s luminosities, as justified in Section [5] It is 
seen that the mass function is built by merging two 
regimes. The highest mass (luminosity) groups are 
seen across the full range of velocities (distances). 
However high mass groups are rare in the nearer 
part of the survey because of the limited volume. 
By contrast, the lowest mass groups are not seen in 
the outer shells but within the range that they can 
be seen they are present in substantial numbers. 

In Fig. Q3] there are overlays of the mass func¬ 
tions built in the two regimes: circles show the 
mass function with contributions from the full 
range 3,000—10,000 km s” 1 and squares show the 



Fig. 14.— Group mass function. Corrections are 
made for loss of information with distance. At 
masses above 3ft” 1 x 10 13 M© the mass func¬ 
tion is built with contributions over the full range 
3,000—10,000 km s” 1 while at lower masses the 
construction is limited to contributions over the 
range 3,000—6,500 km s” 1 . The model fit is 
based on the Sheth-Tormen (ST) formalism with 
9 matter = 0.2 and a count normalization. 


mass function restricted to 3,000—6,500 km s” 1 . 
The numbers for the shallower sample are geomet¬ 
rically augmented to compensate for the reduced 
volume being explored. The two realizations of 
the mass function overlap satisfactorily in the mid 
range of 10 13 — 10 14 Mg. In the plot, circles are 
black and squares are grey above 3ft” 1 x 10 13 M 0 
while circles are grey and squares are black below 
this mass. At the high mass end the contributions 
from the two regimes are consistent but the bin 
VN error bars are larger for the shallower sam¬ 
ple. At the low mass end there are differences 
attributable to completion effects. 


The dashed curve in Fig.Udlillust rates a fh H us¬ 
ing th e Sheth—Tormen modification (ISheth fc Tormen 
1999h of the Press—Schechter formulation of 


the mass spectrum of collap sed h alos f rom an 


initial Gaussian field ( Press fc Schechter 1974ll . 
Fits with re cipes by I Jenkins et al.1 ( 200lh or 


Warren et al. ( 2006h are not significantly different 


from the Sheth—Tormen case. The Sheth—Tormen 
and related descriptions of the mass function are 
motivated by Lambda Cold Dark Matter simula¬ 
tions. 


The model fit minimizes scatter in bins of 
log(Number) over the solar mass range 12.8 to 
14.7 in log units, with all bins in that range given 
equal weight. The fit excludes halos with log mass 
greater than 14.7 because of poor statistics and 
halos below log mass 12.8 because of probable 
systematics. Over the fitted range, the scatter of 
0.07 in log counts is achieved with a low value 
of 9 matter ~ 0.2, assuming a topologically flat 
universe. However, while the shape of the model 
fit is in excellent agreement with the data, the 
count normalization is poor. There will be cosmic 
variance. As a measure of statistical fluctuations, 
there are 7% more 2MRS galaxies in the sample 
volume in the north Galactic hemisphere than in 
the same volume in the south. However, the nom¬ 
inal Sheth—Tormen fit required an augmentation 
in counts by a factor 4.6 to achieve the fit shown 
in Fig. 1141 This unsatisfactory situation requires 
an investigation that goes beyond the scope of this 
paper for reasons discussed next. 

The fit clearly fails at group masses below 
6ft” 1 x 10 12 Mq. There are several possible rea¬ 
sons. There is immediate concern related to the 


4 Kindly provided by Julien Carron and Melody Wolk. 
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luminosity corrections with distance. The adjust¬ 
ment reaches a factor two at 9,000 km s -1 and, in 
the extreme, the adjustment can be as large as 2.6 
or a displacement of 4 bins in log mass. Certainly 
an adjustment is required. Adjustments cause a 
binning uncertainty at the level of ±1 bin. A re¬ 
lated but greater concern is whether the correction 
for missing light should strictly be assigned to one 
halo. The assumption to do so is probably valid 
at high masses but not valid at low masses. High 
mass groups contain a lot of galaxies and statisti¬ 
cally the galaxies missed will be associated. Low 
mass groups contain few galaxies, possibly only 
one that makes the 2MRS 11.75 catalog. It is 
likely that frequently all or part of the added lu¬ 
minosity lies in separate uncatalogued halos. The 
consequence would be that the present procedure 
preserves total light but assigns the light to fewer 
halos than actually exist. There would be deple¬ 
tion of the numbers in the lowest mass bins and 
possible augmentation of numbers in intermediate 
mass bins. 

One wonders about the effect of the loss of low 
surface brightness galaxies in the 2MASS samples. 
The problem is not so much with the systematic 
loss of light. If that light had been captured then 
the coefficients of the correction formula Eq. [7] 
would have been smaller in compensation. The 
greater problem is the loss of halo counts at the 
low mass end where unseen light either at the ex¬ 
tremities of cataloged galaxies or galaxies that are 
entirely missed impacts the mass function. 

It should be possible to evaluate these con¬ 
cerns by ’’observing” mock catalogs. There are 
plans to carry out this experiment. However ul¬ 
timately the 2MRS 11.75 catalog is not optimal 
for a study of the faint end of the mass function 
because of the loss of low surface brightness sys¬ 
tems. As was discussed in the companion paper 
T15, it is apparent that there is an abundance 
of groups with low mass that are entirely missed 
in the 2MRS 11.75 compilation and would con¬ 
tinue to be missed with 2MASS samples to fainter 
limits. Surveys that target low surface bright¬ 
ness sys tems can help clarify the situation at low 
masses (van den Berghlll959: Fisher fc Tullvlll981 


2005j). An anticipated future paper in this series 
will present a group catalog for galaxies within 
3500 km s -1 that samples low surface brightness 
systems. That catalog of nearby groups will give 
much better definition of the low mass end of the 
mass function. 

There can readily be an improvement of the 
observed mass function at the high mass end 
as well. The 2MA SS Extended Source Catalog 
( Jarrett et al. 2000h augmented by Tom Jarrett 
(private communica tion) with velocities from the 
6dF Galaxy Survey ( Springob et al. 2014 ?) and the 
Sloan Digital Sky Survey extends a factor 2 greater 
in depth over 3/4 of the sky, giving access to a fac¬ 
tor 6 more volume. 

While group catalog extensions to the very local 
volume and to larger distances will greatly clarify 
the observational status of the group mass func¬ 
tion, there remains an issue with the normalization 
of the Sheth—Tormen mass function. There has 
to be a clearer understanding of whether masses 
in simulations and the observed universe are com¬ 
mensurate. This issue will be given attention in a 
future publication. 

7. Summary 

John Huchra’s legacy, the 2MASS Redshift Sur¬ 
vey with almost all-sky completeness and photo¬ 
metric integrity, provides the best available de¬ 
scription of the redshift dist ribution of galaxies 


in our corner of the universe (Huchr a et al 
Erdoedu et all[20061 iNuza et ah 20141) . The 


2012 ; 


near- 


infrared photometry captures light from the domi¬ 
nant baryonic component of galaxies, the old stars, 
with minimal loss from obscuration. Variations of 
this sampl e have already been used to build group 


Karachentseva&J£ajachentsev[l998; 


1999 : Karachentseva fe; Karachentsev 2000l) . Wide 


Karachentseva et 


field neutral Hydro gen surveys provide comple- 
mentary inventories ( Mever et al.il2004 : Giovanelli et al 


catal ogs (ICrook et al.1 l2007t lLavaux fe Hudson 
12011 ). In the current paper, in addition to us¬ 
ing the full 2MRS 11.75 catalog, the groups are 
constructed to conform to scaling laws defined by 
detailed studies of individual groups with masses 
10 11 to 10 15 M 0 , discussed in T15. 

While groups have been defined over the full ve¬ 
locity range of the 2MRS 11.75 catalog, the near¬ 
est and furthest groupings are suspect. Due to 
ah the loss of low surface brightness galaxies from 
the 2MRS catalog, this sample is not optimal for 
the construction of low mass groups, best stud¬ 
ied nearby. At high redshifts, once only galaxies 
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brighter than the characteristic magnitude of 
the exponential cutoff are being accessed, the cor¬ 
rection factor for missing light creates unaccept¬ 
able uncertainties. The most useful domain of the 
present catalog is 3,000 — 10,000 km s _1 in the 
CMB frame. 


In this shell from 3,000 to 10,000 km s -1 , the 
present catalog contains 24,044 galaxies in 3,461 
groups of two or more and 10,144 singles. The 
Perseus Cluster, Abell 426, is the most popu¬ 
lated group with 180 assigned members. There 
are 5 groups with masses greater than 10 15 Mq 
(A2199, Coma, Ophiuchus, Norma, and Perseus), 
182 groups in the next decade of mass down to 
10 14 Mq, and 3514 more groups in the subse¬ 
quent decade down to 10 13 Mq. The densest 
concentration of clusters are in the regio n histori¬ 


cally called the ’’Great Attractor” (Dressier et al 


198 7). which is_ th e core of the Laniakea Super¬ 


cluster ( Tullv et al . 20141) , and the Perseus—Pisces 
chain ( Haynes fc Giovanelli 1988T) . 


The construction of the group catalog all the 
way down to ’’groups” of one permits the con¬ 
struction of the mass function of groups. The 
volume considered is between 3,000 and 10,000 
km s _1 in the CMB frame, absent 9% of the 
sky in the zone of obscuration. Every 2MRS 
11.75 galaxy in this volume contributes to halos 
as small as 2 h _1 x 10 11 Mq and as massive as 
2h~ 1 x 1 0 15 Mq. The halo fitting function deter¬ 
mined by Sheth fc Tormenl ( 19991 ) provides a good 
fit to the shape of the observed mass function with 
the choice Q matter ~ 0.2 in a flat universe. How¬ 
ever, the count normalization is poor. This situ¬ 
ation requires further analysis but first there are 
available ways that the observed mass function can 
be improved with catalog extensions to the local 
volume and to greater distances. 


The mass in halos in the 3, 000 —10, 000 km s _1 
shell can be summed. Divided by the volume, the 
product is a density that can be compared with 
the critical density for a closed universe without 
vacuum energy. This normalized density in bound 
halos is ^collapsed = 0.16 ± 0.02. The error is 
the quadrature sums of three components: cosmic 
variance gauged by the north—south difference of 
±0.009, 20% uncertainty in the amplitude of the 
lost light correction factor resulting in a summed 
fractional uncertainty of 0.014, and a 10% uncer¬ 
tainty in the conversion from light to mass, an 


uncertainty in Vtcoiiapsed of 0.016.. 
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